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The contour integrals, occurring in the arbitrary-order phase-integral quantization conditions 
given in a previous paper, are in the first- and third-order approximations expressed in terms of 
complete elliptic integrals in the case that the charges of the Coulomb centres are equal. The eval- 
uation of the integrals is facilitated by the knowledge of quasiclassical dynamics. The resulting 
quantization conditions involving complete elliptic integrals are solved numerically to obtain the en- 
ergy eigenvalues and the separation constants of the lsa and 2pa states of the hydrogen molecule ion 
for various values of the internuclear distance. The accuracy of the formulas obtained is illustrated 
by comparison with available numerically exact results. 
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1. INTRODUCTION 

In a previous paper JlJ the general two-centre Coulomb problem was treated according to the phase- integral method, 
briefly described in the appendix of that paper, and arbitrary-order quantization conditions were given, valid uniformly 
for all energies. 

The symmetric case, where the charge numbers Z\ and Z 2 of the two Coulomb centres are equal, represents for 
Z\ = Z% = 1 the hydrogen molecule ion H+. This case has for natural reasons been studied more extensively than the 
general two-centre problem. Thus, the ion was the subject of treatments already with the aid of the old quantum 
theory [|j and in the early days of quantum mechanics Q . The simplicity of the hydrogen molecule ion grants it an 
analogous unique position in molecular physics as the hydrogen atom possesses in atomic physics. For instance, 
plays a role of fundamental significance in the theory of chemical bonding, and it is also of importance in the study 
of stellar atmospheres. 

For the background of the present paper we refer to [Q. In sec. 2 of the present paper we present the quasiclas- 
sical dynamics of the associated particle system by setting up the equation of motion for a particle moving in the 
relevant potential and express the solution in terms of Jacobian elliptic functions. In sections |^ and ^ we express the 
quantization conditions given in [jl| in terms of complete elliptic integrals by using the solution in sec. 2. We intro- 
duce, analogously as Lakshmanan and Kaliappan ||, Lakshmanan et al. |J and Lakshmanan et al. 0, convenient 
transformations to elliptic functions, in order to be able to express both the real and the complex contour integrals, 
occurring in the quantization conditions, in terms of complete elliptic integrals. In choosing our transformations we 
exploit the symmetry of the functions R(j]) and Q 2 {rf), introduced in jlj], already from the begining of the calculations, 
since this is much simpler than to particularize formulas for the general case of arbitrary nuclear charge numbers Z\ 
and Zi to the case Z\ = Zi. Thus we obtain simple expressions for the quantities appearing in the quantization 
conditions in O. It should be remarked that if one particularizes the ^-quantization conditions involving complete 
elliptic integrals that are valid for arbitrary charge numbers Z\ and Zi (to be derived in the subsequent paper Q|) to 
the case Z\ = Z%, one must in general make further cumbersome transformations in order to bring the quantization 
conditions in question into the same form as those obtained by assuming from the begining that Z\ = Zi. The 
functions -R(£) and Q 2 (£), introduced in [|), do not display an analogous symmetry as R(rj) and Q 2 {rf) when Z\ = Zi. 
The treatment of the ^-equation is thus the same whether Z\ = Zi or Z\ 7^ Zi, and it is analogous to the treatment 
of the ^-equation in the general case when Z\ may be different from Zi, which is treated in the third paper Q in 
our series of papers concerning the phase-integral treatment of the quantal two-centre Coulomb problem. All the 
contour integrals needed in the calculations are expressed in terms of complete elliptic integrals, corresponding to 
the use of the first- and third-order phase-integral approximations. Since complete elliptic integrals can easily be 
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evaluated by means of standard computer programs, computational difficulties that may arise in direct numerical 
calculations (especially of higher-order terms) are thereby eliminated. One achieves also the possibility of being able 
to utilize well-known properties of complete elliptic integrals, for instance series expansions, for analytic studies on the 
basis of the quantization conditions. Some of the first-order contour integrals have been given in terms of complete 
elliptic integrals previously by Strand and Reinhardt 0], but to the best knowledge of the present authors, even in 
the first-order approximation only partial results have been given. 

The phase-integral formulas derived in this paper are quite flexible due to the presence of two arbitrary parameters 
C and C introduced in the base functions Q(rf) and Q(£) in eqs. (3.2b) and (3.2a), respectively, in In accordance 
with the discussion in sec. 3.1 of Q it is appropriate to choose C and C such that the first-order results are exactly 
equal to the third-order results, and we confirm in sec. | the convenience of this approach by investigating for the Isa 
and 2pa states of the hydrogen molecule ion the accuracy of our phase-integral quantization conditions. 



2. QUASICLASSICAL DYNAMICS 



We have seen in sec. || of |jj that the quantization conditions require the evaluation of various contour integrals. 
In the present section we shall restrict ourselves to the ^-equation. The shape of the contours depends on the choice 
of A and on the real or complex nature of the zeros of Q 2 (f]); see Figs. 2 - 4 in [Q. 

Letting, when —Q 2 (i]) is a double-well potential, the_zeros of Q 2 (rf) be 771,772,773,774 when A = \m\ ^ and 772,773 
when A = 0, we shall evaluate the integrals a = f3 and K in the following cases separately. 

1. A = \m\ ^ : 

(a) Subbarrier case (771 = —774 and 772 = —7/3 are real) 

(b) Superbarrier case (771,774 are real and 772,773 are complex conjugate) 



2. A = 0: 

(a) Subbarrier case (only two zeros, 772 — —773; they are real) 

(b) Superbarrier case (only two zeros, 772 and 773; they are complex conjugate) 

In each one of these cases the evaluation of the contour integrals in terms of complete elliptic integrals is facilitated 
by the knowledge of the associated quasiclassical dynamics. In particular the relevant substitution for the integration 
variable in the various contour integrals can be obtained with the help of, for example, Byrd and Friedman |jTo| . 
However, the meaning of these substitutions can be attributed to the associated solution of the equation of motion of 
the corresponding classical problem as in the case of three-dimensional anharmonic oscillators [Q. 

To illustrate what has been said above, we shall consider the subbarrier case of A = \m\ ^ 0. (For the other cases 
the analysis can be carried out in a similar way.) The base function Q{rf) for this case is 



where a — 774 = —771, 6 — 773 = —772 and p 2 = —^Er\ 2 . Let us define a new integration variable u by writing 



u 



(i-v 2 )Q(v) : 



(2.2) 



then 



^ = (1 - 77 2 )Q(77) = Q(77). (2.3) 
au 

Before solving (2.3), we shall relate this differential equation to the equation of motion for the quasiclassical motion 
of a particle with a mass /j in a certain potential by differentiating (2.3) with respect to u, getting 

du 2 2 dr) ' 1 ' ' 
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Defining now the "time" variable u as 

u = iiy/Ji, (2-5) 

we obtain the equation of quasiclassical motion 

& = (^a) 
aw 1 an 

U(V) = ~\q 2 (ii)- (2.6b) 

We can consider (2.6a) as representing the motion of a particle in the potential U{rj). We shall now solve the differential 

equation fl2.3|) in terms of Jacobian elliptic functions when Q 2 (n) is given by (^jj) . The potential U (rj) is then given 
by 

u(v) = -\q 2 (v) = -^-v 2 ) 2 Q 2 (v) 

= -\ P 2 (a 2 -V 2 )(v 2 -b 2 ) 

= lp 2 (v 2 -a 2 )(v 2 -b% (2.7) 



and according to (2.3) and (2.1) 



^= P V(a 2 -V 2 )(v 2 -b 2 ). (2.8) 
du 



By solving (2.8) with respect to n we obtain 

= lPf b _ r (2.9) 
an [f(u - u )\ 

where uq is an integration constant, and the modulus k of the Jacobian elliptic function is given by 

k 2 = l-^. (2.10) 
Similarly we can for all other cases interpret the corresponding quasiclassical dynamics by the above type of treatment. 

3. CASEA=|m|/0 

In order to express the contour integrals, occurring in the quantization conditions pertaining to the ^-equation and 
the ry-equation, in terms of complete elliptic integrals, we transform the integrals in question from the £- or 77-plane to 
another complex plane, the u-plane, on which the Jacobian elliptic functions cnu, smt and dnu are defined. We make 
frequent use of formulas in |^p| . As mentioned above the transformation in question can be attributed to quasiclassical 
dynamics. 

3.1. Four real zeros of Q 2 {n) and Q 2 (£) 
3.1.1. The quantities a — (5 and K pertaining to the n-equation: Subbarrier case [Fig. 3(a) in Ref. 1] 
Putting 774 = —rii = a and 773 = —7/2 = b, we have the base function (2.1), that is, 

QM=p v^w -n\ ,3,, 

1 — 77 z 
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Using the appropriate transformation on p. 54 in E3|, or equivalently the quasiclassical solution given in sec. II, we 
put 



dnu (l~k 2 sn 2 u)' 



(3.2) 



Noting that the loop 772 — > ?7i — > T)2, that is —b — > —a — > —6, in the ?7-plane, denoted by r_;,._ a , represents the contour 
A Q in Fig. 3(a) of (T) which corresponds in the u plane to — » AT — > 2 AT, where K is the complete elliptic integral of 
the first kind, we obtain for the first-order contribution to a the following transformation of the original integral in 
the 77-plane into an integral in the it-plane: 



v(l) - 



pi 



2 A' 



rZ// 



<? 2 7 \du J 1 — T] 2 
k% 2 1 /- 2K sn 2 w(l-sn 2 w) 



1 — 6 2 2 7 (1 — fc 2 sn 2 u)(l — v 2 sn 2 u) ' 
which, after decomposition of the integrand into partial fractions and use of recurrence formulas in [hOl, yields 



(3.3) 



a 



(i) 



E{k) 



1 



K(k) + k 2 1 



1 



n(z/ 2 ,fc) 



where 



a 2 -b 2 fc 2 _ 1 2 

a 2 (l-6 2 ) ~ l-& 2 '- 9 ~ a' 



(3.4a) 



(3.5) 



Here K(k), E(k) and II(^ 2 , k) are complete elliptic integrals of first, second and third kind, respectively. Similarly we 
obtain for the third-order contribution to a (see eqs. (3.13b), (A. 5b), (A. 6b), (A3), (2.9b) and (3.2b) of jy), 



a 



(3) 



C + - 



1 



1 



1 - J] 2 ) 2Q(1 - rj 1 



(dQ(rj)\ 
V dr, J 



dr), 



(3.6) 



where C is the parameter introduced in the base function Q(rj) in eq. (3.2b) of JjJ. After evaluation of the integrals 
we obtain 



a 



(3) 



i^K(k) + 9 

2p { ' 2p(l - b 2 )v 2 



'k 2 K(k) + [y 2 - fc 2 )n(^ 2 ,fc)] 



.9(1 - b 2 ) 

8,pb 2 k 4 



[PxK(k) + P 2 E{k) + P 3 n(v 2 , k)}, 



where 



Pi = i[-9fc 4 + fc 2 (8 + 5^ 2 ) + 4^ 2 - 8], 
P 2 = llk 2 (-4-v 2 ) + (8-^ 2 )} 



(3.4b) 

(3.7a) 
(3.7b) 



and 



2\2 



P 3 = A{u 2 - k 2 ) 



(3.7c) 



Analogous calculations can be performed to evaluate the quantity K. For this purpose we make use of the appro- 
priate transformation on p. 58 in |10|, that is 



2 1,2 2 
rj = b sn u. 



(3.8) 
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The first-order (see eqs.(3.15b), (A5b) and (A6a)_of 0) and the third-order (see eqs. (3.15b), (A5b), (A6b), A(3), 
(2.9b) and (3.2b) of @) contributions to K (= irK) are 



and 



where 



i 
2 

pis* 



Q{r])dr] 



a k 

2 i-IK 



g Jo 1 — v 1 sn 2 u 



-du 



E{k) + fr 1 ( 1 - — 



iT(fe) + K - * ) 1 



noAfe) 



ttKi = - 



g (3) (77)^77 



-21m- 



-C 



1 1 n -3c 
1 - r; 2 ; 2Q(1 - 7 ? 2 ) ~ 8° (??) 



d?7 



%r(*)-2n(^*)+ » 

p p Ab^p 



1 



3A- 



^ (-3^ 2 fc 6 - fc 4 + 8^ 2 fc 2 - 7z^ 2 + l)E(k) 



v 2 = b 2 



9=-, k 2 = b - 2 , k 2 = l-k 2 = 



(3.9a) 



The integrals a' and /3' for the contours A a > and A^/ in Fig. 3(a) in |Q are obtained from the formulas a' 
and /?' = /?+ ^; see eq. (3.18a) in §. 



(3.9b) 
(3.10) 

„ i Air 



3.1.2. The quantities L and L' pertaining to the (^-equation [Fig. 1(a) in Ref. 1] 



Denoting the four real zeros £i < £2 < 1 < £3 < £4 of Q 2 (£) by the simpler notations d<c<l<b<a, respectively, 
used in pi, we have 



Q(0=P 



[( a -fl(g-6)(g-c)(g-d)]i 



Using the appropriate transformation on p. 120 in M], we obtain (cf. sec. II) 



- b 



1 — isfsn 2 u 
a 



< 1. 



(3.11) 



(3.12) 



Noting that the loop £3 — > £4 — > ^3 , that is b — > a — > 6, in the £-plane, denoted by I\ a , represent s the contour A^ 
in Fig. 1(a) of jij and corresponds to — > if — > 2if in the w-plane, and using the transformation (3.12), we obtain 
the first-order contribution to L through the following transformation of the original integral in the £-plane to the 
it-plane: 



jl f 2K (d£) 2 du 

" 2~g] \du) ~e^l 

2p{v 2 -v 2 )(y 2 ~v 2 1 ) f 2K sn 2 u{\- sn 2 u){\ - k 2 sn 2 u) 

g Jq (1 — v 2 sn 2 u) 2 {l — v 2 sn 2 u){\ — v 2 sn 2 u) 
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where 



g = 



[(a-c)(b-d)Y- 



(a-b)(c-d) 
(a-c)(b-d)' 



(3.15) 



Note that a, 6, c, d, and hence also V\, ^2, ^3, 5, fc, depend on the choice of the parameter C in the base function Q(£); 
cf. (3.2a) in Decomposing the integrand in (3.13) into partial fractions, and using recurrence formulas in we 
obtain the final formula 



2p 



3k 2 

1 _ 2fc 2 + -j- ) K(k) - 3E(k) 



3k< 



+ I 2(1 + fc 2 ) - v\ - 11(1/?, k) - QSi 
where we have introduced the "universal" function and 



i=l 



Ci = 


2[2v 2 v 2 -v 2 V 2 


-44] 


(i/f - j/ 2 )(i/ 2 


-4) 


c 2 = 


2^f-^ 2 




[y\ - v\)(y 2 - 


4Y 


c 3 = 


24(4 - 4 




(4~4)(4~ 


4Y 



(3.16a) 
(3.16b) 

(3.17a) 
(3.17b) 
(3.17c) 



1 



[(4 + 24k 2 - 3k 2 )K{k) - (4 + 4k 2 - 3k 2 )E{k) 

1 = 1,2,3, 



3fc 2 
3k 2 

+-r(i-4)(k 2 -4)m4,k)-K(k)} 



which can also be written as 



Si = T 



,K(k)-E(k) (3-2v 2 )E(k) 



k 2 



k' 2 



311 



i = 1,2,3, 



the last formula being valid if v 2 and k 2 fulfil the conditions stated in section 117.03 in [ p~0| 
Similarly we get for the third-order contribution to L: 



c 



1 



1 



e-ij 2Q(0(e-i) * 1^ 



that is 



Z< 3 > =-HW(v 1 , V2 ,v 3 ,g,k,C) 



GM4 - 


4)(4 


-4) 


9_ 


C + 


4(4 - 


-4) 


Ap 


±4(4 


-4) 



Y,ClK(k) + J2 D i E (k) 



i=i i=i, 

,,2/„2 ,,2^, 



,2/, ,2 



K(k), 



(3.18) 



(3.18') 



(3.16c) 



(3.16d) 
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where is another "universal" function, C is the parameter in the base function Q(£) (cf. eq. (3.2a) in |Q), and 
4 



3k' 2 
4 



jj [fc 4 + 2fc 2 - 2 + (2^ + v\ + i/f )(1 - 2k 2 ) + {2v\(v 2 + v 2 ) + i/fz/f + z/ 4 }(4 - 3fc 2 )] 

11 - 10fc 2 



3fc' 2 



fc 2 



fc 4 



(14- 6fc 2 - 7fc 4 ) 



C 2 = 


4 


2 2 
^3 


c» = 




3fc 4 




16 - 



4 4 . 4 4. 4-4- p. 4 v v 

2^ 2 ^ 3 +^1/3 + I^IA. - 2^^^ 



4,, 4 



,4,, 2,, 2 



(3fc 4 - 6z^fc 2 + 2v\ + k l v\), 



2v 2 v 4 v 4 ( v\\ 4v 4 v 4 v 4 



k 2 



1^1 

2 -y 



3fc 4 



,2,2 , ,,2, ,2 



1 2 3 / o i i 2 \ i 1 / 2 2 2 2 i 2 2\ 

fc4 ( 2 + fc ) + -piTi v i +4^2^) 



(3.19a) 

(3.19b) 

(3.19c) 
(3.19d) 



D 1 = — [2- 3fc 2 - 3fc 4 + 2fc 6 - {2vl + v 2 + j/f)(l - 4k 2 + fc 4 )] 



D 



3k' 4 
4 



3k 
4 



3fc 



(11 - 20fc 2 + llfc 4 ){2^ 2 M + vVA + "3)} 
k 2 



fc 4 



-(14- 13fc 2 - 13fc 4 + 14fc 6 ) 



^ 2 - 3fc 4 
D 3 = D 2 , 



16 ^ 2 [3fc 2 -, 2 (l + fc 2 )], 



k 2 



0222 

,-,,-> 1 ,,,2, .2, 1 1 2 3 / -J 1 ^,2-j 



"K^2 f+ 4^2^.3) + 



fc 2 



-(1 



(3.20a) 

(3.20b) 
(3.20c) 

(3.20d) 



The integral L' for the contour A^, in Fig. 1(a) in Q is obtained from the formula L' = L + ^-n. Therefore L W 
and L ^ can be obtained from (3.16a,b) and (3.16c,d). 



3.2. Two real and two complex conjugate zeros of Q 2 {ri) and Q 2 (£) 
3.2.1. The quantities a = f3, K, L and L' pertaining to the n-equation: Superbarrier case [Fig. 4( a ) ° r Fig. 2 in Ref. 1] 
Putting 7/4 = — 771 = a, i]2 — r/3 — —iai, we have 

[(a 2 - r/ 2 )(a 2 + rj 2 )]^ 



Q(v)=p- 1 

1 — r/ z 

Using the appropriate transformation on p. 133 in pi, that is, 

,2 



2 2 

a an u, 



(3.21) 



(3.22) 



we now utilize the fact that the Jacobian elliptic functions are doubly periodic, one of the periods being complex. 
Thus, the loop 771 — > 772 — * 771, that is —a — > —ia\ — > —a, in the ?7-plane, denoted by r_ a ._i Ql , represents the contour 
A a in Fig. 4(a) of and it corresponds in the w-plane to — > A" + iAT' — » 2A~ + 2iAT' where A'' = K(k'). Denoting 
by 2G( 2n+1 ) the integral occurring m the definitions (3.13b) and (3.156') in § when Z\ = Zo. we have the following 
first-order expression: 



2 JA„ 
1 

2 



Q(v)dr) 

Q(v)drj, 



(3.23) 
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which when transformed to the variable u becomes 



G^ 1 



2 r 2K+2iK' 



pa 

~2g J Q 



sn udn u 
1 — a 2 cn 2 u 



-du. 



After evaluation of the integral in ( 3.24 ) we obtain 



where 
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E(k) + ^(1 - a 2 )K{k) - (% + k' 2 ) Il(v 2 ,k) 
-h, \E(k') + (k 2 + a 2 k' 2 ) [n(l - v 2 , k') - K(k')] } , 



(3.24) 



(3.25a) 



Similarly one obtains 



Q(V=-^l K (k) 
2p 



8v 2 p 



a 2 - 1 



, 9 = 



(a 2 + a 2 ) 2 



- k 2 - 
-, re — 



a 2 + a 2 



, k 2 = 1 - k 2 



2p(l - a 2 ) 



(y 2 - i)g 

4p 



n(v 2 ,k)-*-l a2{1 - a2) 



2 V k' 2 a 2 + k 2 



(-4fc 2 + Av 2 - l)K(k) + —(-8k 4 + 8fc 2 (l - v 2 ) + 1 + 7v 2 )E(k) 

3k 1 



2p(l - a 2 ) 



(^ ~ l)g 
4p 



(1 - a 2 )K(k') - i/ 2 n {v 2 - 1, fc') - - 
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8v p 



\{-Ak 2 + Av 2 - l)K(k') + -I_ (-8fc 4 + 8fc 2 (l - ^ 2 ) + 1 + 7^ 2 )[K(fc') - 



(3.26) 



(3.25b) 



Recalling the definition of G' 2n+1 - ) above (3.23) in the present paper, one sees that according to (3.13b) in p| the 
first- and third-order contributions to a are 



= ReG« 



(3.27a) 



«(3) =ReG( 3 ), 

and that according to (3.156') in Q the first- and third-order contributions to if (= tt-K") are 

ttK = -2ImG (1) , 



(3.27b) 



(3.28a) 



irK 2 = -2ImG (3) . (3.28b) 

The integrals a' = f3' associated with the contours A a / and A^/ in Fig. 4 in Q are obtained from the integrals 
a = /3 by means of the relations (3.18a) in 0, that is, a' = /3' = a + 4? = /3 + 4p. The integrals 1/ and L' associated 
with the contours A^ and A^/ in Fig. 4(a) in [|l] can be obtained from the integrals a = j3 by means of the formulas 
L = a + [3 = 2a and V = L + \m\. 

The only essential difference between —Q 2 {rj) in Fig. 4(a) in Q and — Q 2 (ry) in Fig. 2 in H is that in the former 
figure there is an underdense barrier of —Q 2 {rj), while in the latter figure —Q 2 (ji) has a single minimum. For the 
case in Fig. 2 one has therefore the formula L { - 2n+1 ^ = a( 2n+1 ) + /?( 2 ™+ 1 ) = 2a( 2 ™ +1 ), with the expansions (3.27a) 
and (3.27b) for a^ and a^ 3 \ and the formula L ( 2 ™ +1 ) = L,i' 2n + 1 ) -f |m|£ m ,o- The case in Fig. 2 has, however not 
appeared in our applications. 
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3.2.2. The quantities L and V pertaining to the ^-equation [Fig. 1(a) in Ref. 1] 

Denoting by £3 = b and £4 = the real zeros of Q 2 (£), and by c and c* the complex conjugate zeros £1 and £2 of 
Q 2 (£), we have 

Q(0=p j ■ I 3 - 29 ) 

Defining 

c = b\ — iai, c* = bi+iai, (3.30) 



A = [(a-b 1 ) 2 + a 2 1 ]i, B=[(b-b 1 ) 2 + a 2 }^ (3.31) 

and using the appropriate transformation on p. 133 in (To), we get 

aB + bA + (bA - aB)cnu 
«- A + B + (A-B)cnu ' ( ' 

Noting that the loop b — > a — > b in the £-plane (denoted by I\ a for the contour A^) corresponds to the path 



— > 2_ftT — > in the u-plane, and using the transformation (3.32), we obtain for the first-order contribution to L 
the formula (see Fig. 1(a) in [Q), 

JA L 

= Pfr-'*)^-'*) t K snWu 

9 Jo (1 + v\cnuy(\ + vicnu)[\ + v^cnu) 

where 



_ A-B _ (l + b)A-(l + a)B _ (l-b)A-(l-a)B 
Vl ~A + B> V2 ~ (l + b)A+(l + a)B> y *~ (l-b)A+(l-a)B> {3M > 

- ^ {a - b? 7i A - B? - (3-35) 



r AB ^AB 



By evaluating the last integral in (3.33), and introducing a new "universal" function H^ l \ given by eqs. (2.19), 



(2.20a-c) and (2.21) in [§ with j = 0, we obtain 

L (1) = -2ReiJ (1) (^i,^ 2 ,^3,5,fc,C')[with j=0] (3.36a) 



2£ J-L 

(^ 2 - 2k 2 )K(k) - 2v\E(k) + v x kix\ + CiJi \ , (3.36b) 



where 
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Co = 



(v X - V 2 )(V2 - "3)' 



(3.37b) 



C, 



(l/l -V2)l4 



{V\ - ^)(^3 - VnY 



(3.37c) 



Jj = Re5j(with j = 0) 
1 



fc^i - ^)x(fc) + vfE{k) - {e + v 2 k' 2 )n ( -fi— , fe 



vt-V 



2kv? 



7r,i = 1,2,3; 



cf. for the definition of 5j eq. (2.21) in |q]. The third-order contribution to L is 



L (3) = -2ReF (3) (^i,j/ 2 ,^3,3,fc,C')[with j=0] 

^ 2 (^2 - V3) 2 
^2^3(^1 - ^3)(i / 2 - Vl) 



(3.38) 



(3.39a) 



— ^ (4C- 



K(k) 



[XK{k)+YE{k)\ 



(yx - v$)(vx - v 2 ) 

where is another "universal" function, given by eqs. (2.22) and (2.23a-c) in || with j — 0, and where thus 



(3.39b) 



X = - 



(l + 4fc 2 ) (3 + 4/fc' 2 ), 2 , I' 2 , 2l2 

7, 1 n \ V l + 2l/ 1^2 + 2^1 1S 3 + - X77(17 - 4fc )^x^2^3 

3 3 3fc z 



y 



2vi(vi +V 2 + V3) + 2^2^3 + v\ f — + — ] , 

V ^3 ^2 / 

= ^ (1 + 8fc 2 - 8fc 4 ) + ^(2fc 2 - 1)K 2 + 1vxv 2 + 2v x v z + v 2 v z ) 



3fc 2 



(17- 



(3.40a) 



(3.40b) 



The integral V associated with the contour A^, in Fig. 1(a) in Q is obtained from the formula V = L + ^-tt. 
One has therefore the formulas Z'W = ZW + ^tt and = with and given by (3.39a) and (3.39b). 

4. CASE A = 
4.1. Two real zeros of Q 2 (n) and Q 2 (0 
4-1.1. The quantities a = (3 and K pertaining to the rj-equation: Subbarrier case [Fig. 3(b) in Ref. 1] 
Putting 7/3 = — r)2 = b as before, we have 



Q(v)=p 



(ry 2 - b 2 ) 



1 — I] 2 

Using (3.2) with a = 1, the first- and third-order contributions to a become 



(4.1) 



= p[E{k)-(l-k 2 )K{k)}, 



(4.2a) 
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v(3) = 



1 



2pk 2 



(1 - Ck 2 )K(k) - E(k) - ( (3k 2 - 8)K(k) - (7k 2 - 



k A 



where k 2 = 1 — b 2 . 

The first- and third-order contributions to K(= irK) are 



nK = 2p[E{k) - (1 - k 2 )K(k)}, 



where k 2 = b 2 . 



(4.2b) 



(4.3a) 
(4.3b) 



4-1-2. The quantities L and L' pertaining to the ^-equation [Fig. 1 in Ref. 1} 
If the real zeros of Q 2 (C) are £3 = c (< a ) an d £4 = a i we have 

Q(0=pf K " C)( °~ l* (4 4) 

We shall treat the three cases 1 < c < a, — 1 < c < 1 < a and c<— 1 < 1 < a separately. 

Case 1 < c < a [Fig. 1(a) in Ref. 1] 

For this case we use in |L0| the transformation on p. 120 and the formula in section 256.19 with a suitable choice 
of parameters to obtain the first- and third-order contributions to L as 



= J Q(£)d£ 

(a - c)(c- l)g 
P 2v 2 (v 2 -k 2 ) 

m = CgK(k) , g 



2p 



[(2v 2 -v 4 - k 2 )Il(v 2 ,k) - u 2 E(k) - (v 2 - k 2 )K(k)] , 
2t^W P [^ (fc2 " V2)K{k) + W (fc2(fc2 ~ VA ~ ^ + 2vi) m 



where 



—^4 [(-2 + k 2 + 2v 2 - v 4 )K(k) + 2 (k 4 - (1 + v 2 + v 4 )k 2 -v 2 + 2v 4 + 1) E(k)} , 

24v^p(c — 1) 

2 a-c 2 2 2(a-c) 



a-V 



(a-l)(c+l) 



(4.5a) 



(4.5b) 



(4.6) 



[(a-l)(c+l)]» 

To obtain L' one can use the formula L' ( - 2n+1 '> = i( 2 "+ 1 ) + ^5 mfi with LW and L^ 3 ) given by (4.5a) and (4.5b). 
Case -1 < c < 1 < a [Fig. 1(b) in Ref. 1] 

We use in [jl0| the appropriate transformation on p. 120 and the formula in section 256.17 (with b = 1 and d = —1) 
to obtain the first- and third-order contributions to L' (see Fig. 1(b) in jlj) as 



f.'d) = BP 
5 



(4.7a) 



~, (3) = gCif(fc) | ff 

2p Ap>v 2 (\ - c) 



1 - ^ ) + (k 2 (2k 2 - 1 - 2v 2 ) + v 4 ) 



E(k) 



k 2 k' 2 



9 



16(1 — c)pv 2 



— ((2 + v 4 + 2v 2 )k 2 - 3fc 4 - 2v 4 ) K(k) 



3k' 



- ( -(^ + y 1 - i) + fcV - y 1 + 2) + ^ ) 



(4.7b) 
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where 



a- 1 



2 _ (l + c)(o-l) 



2(a - c) 



(4.8) 



Case c <-!<!< a [Fig. 1(b) in Ref. 1] 

Using in |lO| the transformation on p. 120 and the formula in section 256.20 with a suitable choice of parameters, 
we obtain the first- and third-order contributions to L' as 



with 



: '95 



(i) _ {a - c){l - c)gp 



£'(3) 



2u 4 

CgK(k) , g 



p 4pv 2 
fJ 



[-v 2 E{k) + (y 2 + k 2 )K{k) + (i/ 4 - k 2 )U(u 2 , k)] , 



k 2 ' 



2,2pv 2 



+ 



3k 2 k' 4 



{w^[-^ + k\2-^-2 V ^)-^] K{k) 
k 6 + k 4 {v 2 ~ v A ) + k 2 (v 4 + v 2 - 1) - i/ 4 ] 



a- 1 



-j 5 



o+l [(a + l)(l-c)]! 



( g -l)(- c -l) 
(a + l)(l-c) ' 



(4.9a) 



(4.9b) 



(4.10) 



4.2. Two complex conjugate zeros of Q 2 (i~i) 
The case of two complex conjugate zeros of the square of the base function occurs only for the 77-equation. 

4-2.1. The quantities a = j3 and K pertaining to the n-equation: Superbarrier case [Fig. 4(b) in Ref. 1] 
With 7]2 — — ioi and 773 = ia\ we have 



1 — rj z 



a\+rj 2 \ 2 

1 — T) 2 



(4.11) 



Specializing to the case a = 1 in (3.27a,b) and (3.28a,b) along with (3.25a,b), we obtain the first- and third-order 
contributions to a as 



and the first- and third-order contributions to K(= ttK) as 

nK = -2ImG (1) , ttK 2 = -2ImG (3) , 



where 



(5(3) = 9_ 
2p 



G« = ?-{E{k)+i[K{k')-E(k')}}, 
9 



(1 - C)K(k) - E(k) 1 UkQz) + {Sk2 k , 2 7) E(k) 



+1 — 
2p 



with 



-CK(k') + E(k') - i (W(fc') + (8fc [, 2 ?) {^(fcO - £(*')} 



k 2 



(4.12) 

(4.13) 
(4.14a) 



(4.14b) 



(4.15) 



(l + o?)*' " l + o?' 

The integral L' associated with the contour A^/ in Fig. 4(b) in pi] is obtained from the formula L' = a + (3 = 2a 
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5. ACCURACY OF THE PHASE-INTEGRAL QUANTIZATION CONDITIONS FOR THE lsa AND 2pa 

STATES OF THE HYDROGEN MOLECULE ION 



For the lsa and 2pa states of the hydrogen molecule ion one has to put A = 0. The quantization conditions in 
for the lsa state are (3.5a) with 5 = [Fig. 1(a) in Ref. 1] and (3.9) with s = m = [Fig. 4(b) in Ref. 1] 
when ri2 is sufficiently small, but (3.5b) with s = m = [Fig. 1(b) in Ref. 1] and (3.25b) with s a — sp = m = 
[Fig. 3(b) in Ref. 1] when 1*12 is sufficiently large. The quantization conditions in jjj for the 2pa state are (3.5a) 
with 5 = [Fig. 1(a) in Ref. 1] and (3.9) with s = 1 and m = [Fig. 4(b) in Ref. 1] when rj.2 is sufficiently 
small, but (3.5b) with s = m = [Fig. 1(b) in Ref. 1] and (3.25a) with s a = sp = m = [Fig. 3(b) in Ref. 1] 
when r\2 is sufficiently large. After having expressed these quantization conditions in the first and third order of the 
phase-integral approximation in terms of complete elliptic integrals, as described in the previous sections, we have 
use d these quantization conditions to calculate the energy E and the reduced separation constant A' . In subsection 



5.1 we determine C and C as functions of r%2 such that the first- and third-order quantization conditions give the 
same values of both E and A'. For the values of C and C thus obtained, the choice of the base functions Q{rj) and 
Q(£) is optimal in the sense that the most accurate first-order values of E and A' are obtained, since the first- and 
third-order approximations give the same values of E and A'. In this connection we remark that there are quantal 
systems for which one can obtain exact values of the energy by choosing the base function such that the first- and 



third-order results coincide; see p. 1826 in 11 If] and p. 16 in |12j. In subsection 5.2 we determine C and C such 



that the phase-integral quantization conditions give the numerically exact values of E and A' obtained by Murai and 
Takatsu 13, 14|] and establish that the values of C and C thus obtained are in qualitative agreement with the values 



of C and C determined in subsection 5.1 



5.1. Determination of C(r\2) and G(r\2) such that the first- and third-order quantization conditions give the 

same results 

By determining C and C for each value of ri2 such that the first- and third-order quantization conditions give the 
same value of E as well as of A', we have obtained the results in Table I for the lsa state and in Table II for the 2pa 
state of the hydrogen molecule ion; see also Figs. 1 and 2. The phase-integral values of E and A' in these tables are 
in reasonable agreement with the numerically exact values obtained by Murai and Takatsu [jl3| as is best seen 
from Fig.l for the state lsa and from Fig. 2 for the state 2pa. 

It is seen that in Figs. 1 and 2 there are sometimes some irregularities in the values of C , C, \E — Emt\ and 
\A' — A' MT \ for low ri2 values. It should also be noted that C and C approach the correct limiting value \ as ri2 — > 0. 



5.2. Determination of C(rv£) and Cirvz) such that the first-order phase-integral quantization conditions 

reproduce numerically exact values of E and A' 

By determining C and C for each value of r±2 such that the first-order quantization conditions reproduce the 
numerically exact values of E and A' calculated by Murai and Takatsu |l3|, [l4| we have obtained the values of C 
and C presented in Table III for the lsa state and in Table IV for the 2pa state. In Figs. 3 and 4, we have shown 
the dependence of C and C on ri2 for the states lsa and 2pa, respectively. One can carry out a similar calculation 
using the third-order phase-integral quantization condition also and determine the appropriate C and C. However, 
the numerical analysis becomes too laborious and time consuming and so we have not presented the results here. 

To obtain the numerical results in subsections V.A and V.B a general FORTRAN computer program using very 
rapid library routines was written at the Centre for Nonlinear Dynamics, Department of Physics, Bharathidasan 
University, Tiruchirapalli, India. We have carried out the numerical calculations by Silicon Graphics Power Indigo 2 
XZ Graphics Workstation (R8000, 64bit processor) using FORTRAN 77 compiler. 

Some years ago, a direct numerical integration of the contour integrals in the phase-integral quantization conditions 
for the hydrogen molecule ion was carried out by Fil. lie. Anders Hdkback at the Department of Theoretical Physics, 
University of Uppsala, Sweden. By means of this numerical material it was possible to make valuable checks of the 
correctness of the phase-integral quantization conditions expressed in terms of complete elliptic integrals. 
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TABLE I: For the state lsa of H% the values of C and C have been obtained from the requirement that the first- and 
third-order phase- integral results coincide for E as well as for A' . With the use of these values of C and C, the values of E 
and A' have then been obtained from the quantization conditions that are appropriate depending on whether ri2 is sufficiently 
small or sufficiently large. The numerically exact values obtained by Murai and Takatsu [13,14] are given in the columns called 
Emt and A' MT . 



r\i 


C C 


E Emt E — Emt 


A 1 A 1 A 1 A' 

A A MT A A MT 


Sufficiently small m 


0.6 


0.4091767180 0.4953643560 


-1.618424439 -1.6714847145 0.053060275 


-0.1992407489 -0.1992300000 -0.000010748 


0.8 


0.4431528990 0.4884579844 


-1.545199993 -1.5544800915 0.009280098 


-0.3249142019 -0.327900000 0.002985799 


1.0 


0.4681771744 0.4902081420 


-1.455314042 -1.4517863130 -0.003527729 


-0.4706879721 -0.4759469161 0.005258944 


Sufficiently large ri2 


2.0 


0.5255080000 0.5057903368 


-1.109450173 -1.1026342150 -0.006815958 


-1.385866329 -1.393538844 0.007672515 


3.0 


0.5347955900 0.5146860140 


-0.9146268870 -0.9108961974 -0.00373069 


-2.451736492 -2.458030452 0.00629396 


4.0 


0.5255436520 0.5191506880 


-0.7981985544 -0.7960848837 -0.002113671 


-3.564275870 -3.569090310 0.00481444 


5.0 


0.5095344153 0.5210276591 


-0.7256775451 -0.7244202952 -0.00125725 


-4.673981394 -4.677559936 0.003578542 


6.0 


0.4942331566 0.5213617485 


-0.6794162838 -0.6786357151 -0.000780568 


-5.759202308 -5.761839130 0.002636822 


7.0 


0.4829917740 0.5208505340 


-0.6489619089 -0.6484511471 -0.000510761 


-6.817236218 -6.819239945 0.002003727 


8.0 


0.4758245700 0.5199462800 


-0.6279265372 -0.6275703886 -0.000356149 


-7.854452399 -7.856077820 0.001625421 


9.0 


0.4715076730 0.5189057600 


-0.6125711074 -0.6123065640 -0.000264543 


-8.878337258 -8.879752233 0.001414975 


10.0 


0.4688953184 0.5178568616 


-0.6007859807 -0.6005787289 -0.000207251 


-9.894343735 -9.895643269 0.001299534 


15.0 


0.4641710460 0.5135686450 


-0.5668106207 -0.5667156052 -0.000095015 


-14.93203345 -14.93315205 0.00111865 


20.0 


0.4625834317 0.5107892808 


-0.5500740452 -0.5500142593 -0.000059786 


-19.94891764 -19.94996067 0.00104303 


25.0 


0.4617350800 0.5089200200 


-0.5400488855 -0.5400058008 -0.000043085 


-24.95899430 -24.95998443 0.00099013 
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TABLE II: For the state 2pa of the values of C and C have been obtained from the requirement that the first- and 
third-order phase- integral results coincide for E as well as for A' . With the use of these values of C and C, the values of E 
and A' have then been obtained from the quantization conditions that are appropriate depending on whether ri2 is sufficiently 
small or sufficiently large. The numerically exact values obtained by Murai and Takatsu [13,14] are given in the columns called 
Emt and A' MT . 



ri2 


c c 


E Emt E — Emt 


4' A' A' A' 
/i Si-MT A A MT 


Sufficiently small rvi 


0.2 


0.2633959780 0.246516819 


-0.5020720303 -0.5026800000 -0.00060797 


-2.005074078 -2.004020000 0.001054078 


0.4 


0.2688550280 0.234829974 


-0.5091302928 -0.5107900000 -0.001659708 


-2.018393631 -2.016330000 0.002063631 


0.6 


0.2750685470 0.209525191 


-0.5231745581 -0.5243050000 -0.001130442 


-2.041072959 -2.037690000 0.003382959 


0.8 


0.2826340310 0.1470935880 


-0.5524257935 -0.5427400000 0.009685793 


-2.075482537 -2.069270000 0.006212537 


2.0 


0.3435728404 0.5524319489 


-0.6571837692 -0.6675343922 -0.010350623 


-2.523233873 -2.521958177 0.001275696 


3.0 


0.3915235980 0.5260048459 


-0.7015678961 -0.7014183334 0.000149563 


-3.202554526 -3.196382289 0.006172237 


Sufficiently large ri2 


4.0 


0.4279941950 0.5249345450 


-0.6960011808 -0.6955506394 0.000450542 


-4.029227017 -4.025940635 0.003286382 


5.0 


0.4481616320 0.5238924640 


-0.6777410414 -0.6772916132 0.000449428 


-4.942443752 -4.941274459 0.001169293 


6.0 


0.4591233730 0.5226760960 


-0.6577104664 -0.6573105590 0.000399907 


-5.903631412 -5.903659889 -0.000028477 


7.0 


0.4644579210 0.5214162480 


-0.6394674183 -0.6391288554 0.000338563 


-6.890343283 -6.890997919 -0.000654636 


8.0 


0.4666144860 0.5201789140 


-0.6238870346 -0.6236060156 0.000280969 


-7.889744744 -7.890707161 -0.000962417 


9.0 


0.4671264910 0.5189985500 


-0.6108876012 -0.6106549406 0.00023266 


-8.894778515 -8.895880333 -0.001101818 


10.0 


0.4668770890 0.5178931230 


-0.6000950569 -0.5999010686 0.000193989 


-9.901798495 -9.902954530 -0.001156035 


15.0 


0.4641395400 0.5135689400 


-0.5668036046 -0.5667087290 0.000094875 


-14.93214472 -14.93326111 -0.00111639 


20.0 


0.4625830500 0.5107892800 


-0.5500739820 -0.5500141977 0.000059785 


-19.94891896 -19.94996191 -0.00104295 


25.0 


0.4617350700 0.5089200100 


-0.5400488850 -0.5400058003 0.000043085 


-24.95899432 -24.95998445 -0.00099013 
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TABLE III: For the state lsa of H% the values of C and C have been obtained from the requirement that the first-order 
phase-integral results, obtained from quantization conditions that are appropriate depending on whether ri2 is sufficiently small 
or sufficiently large, coincide for E as well as for A' with the numerically exact results obtained by Murai and Takatsu [13,14] 
and quoted in this table as Emt and A' MT . 

ri2 Emt A' mt C C 



Sufficiently small ri2 

0.2 -1.9286203017 -0.0256900000 0.3148379823 0.3506512286 

0.4 -1.8007540595 -0.0957300000 0.3477075009 0.4296524070 

0.6 -1.6714847145 -0.1992300000 0.3681917193 0.4607730539 

0.8 -1.5544800915 -0.3279000000 0.3807005347 0.4777532944 

1.0 -1.4517863130 -0.4759469161 0.1120058693 0.4885447845 



Sufficiently large m 

2.0 -1.1026342150 -1.393538844 0.5095379703 0.5120864521 

3.0 -0.9108961974 -2.458030452 0.5196342600 0.5202577750 

4.0 -0.7960848837 -3.569090310 0.5131695740 0.5235297064 

5.0 -0.7244202952 -4.677559936 0.5001127420 0.5243749920 

6.0 -0.6786357151 -5.761839130 0.4872379220 0.5239023990 

7.0 -0.6484511471 -6.819239945 0.4776859510 0.5227893710 

8.0 -0.6275703886 -7.856077820 0.4715798270 0.5214462830 

9.0 -0.6123065640 -8.879752233 0.4679041020 0.5200867890 

10.0 -0.6005787289 -9.895643269 0.4656852040 0.5188037440 

15.0 -0.5667156052 -14.93315205 0.4617106320 0.5139609640 

20.0 -0.5500142593 -19.94996067 0.4603989670 0.5109975070 



25.0 -0.5400058008 -24.95998443 0.4597077000 0.5090477000 
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TABLE IV: For the state 2pa of H% the values of C and C have been obtained from the requirement that the first-order 
phase-integral results, obtained from quantization conditions that are appropriate depending on whether ri2 is sufficiently small 
or sufficiently large, coincide for E as well as for A' with the numerically exact results obtained by Murai and Takatsu [13,14] 
and quoted in this table as Emt and A' MT . 

T\i Emt A-mt 

C C 



Sufficiently small ri2 


0.2 


-0.5026800000 


-2.004020000 


0.264492583 


0.2439809316 


0.4 


-0.5107900000 


-2.016330000 


0.2712715280 


0.2275797176 


0.6 


-0.5243050000 


-2.037690000 


0.2795333746 


0.2071325406 


0.8 


-0.5427400000 


-2.069270000 


0.2897524086 


0.1952877316 


1.0 


-0.5648136251 


-2.112417232 


0.3022903142 


0.2116506664 


2.0 


-0.6675343922 


-2.521958177 


0.4118503542 


0.5180809640 



3.0 -0.7014183334 -3.196382289 -0.7005640546 0.5327867225 



Sufficiently large ri2 



4.0 


-0.6955506394 


-4.025940635 


0.4300576700 


0.5304281261 


5.0 


-0.6772916132 


-4.941274459 


0.4476172910 


0.5276798434 


6.0 


-0.6573105590 


-5.903659889 


0.4571245682 


0.5253680730 


7.0 


-0.6391288554 


-6.890997919 


0.4617645750 


0.5234038890 


8.0 


-0.6236060156 


-7.890707161 


0.4636608160 


0.5216941120 


9.0 
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FIG. 1: Plots for the lsa state of the ion of (a) C versus ri2, (b) C versus ri2, (c) \E — Emt\ versus 7-12 and (d) A' — A' MT 
versus ri2, when C and C are determined as functions of r r i from the requirement that the first-order phase-integral results 
coincide with the third- order results. Here E and A 1 are the phase- integral values obtained in Table I, while Emt and A' MT 
are the corresponding numerically exact values obtained by Murai and Takatsu [13,14] and quoted in the same table. There is 
a break in each curve between the regions where the quantization conditions for sufficiently small and sufficiently large values 
of ri2 have been used. 
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FIG. 2: Plots for the 2pa state of the ion H 2 of (a) C versus r 12 , (b) C versus r i2 , (c) |S — £mt| versus ri 2 and (d) A' — A' MT 
versus r i2 , when C and C are determined as functions of m from the requirement that the first-order phase-integral results 
coincide with the third- order results. Here E and A' are the phase-integral values obtained in Table II, while Emt and A' MT 
are the corresponding numerically exact values obtained by Murai and Takatsu [13,14] and quoted in the same table. There is 
a break in each curve between the regions where the quantization conditions for sufficiently small and sufficiently large values 
of ri2 have been used. 
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FIG. 3: Plots for the lsa state of the ion of (a) C versus m and (b) C versus ri2, when C and C are determined as 
functions of ri2 from the requirement that the first-order phase-integral results coincide with the numerically exact results 
obtained by Murai and Takatsu [13,14]. There is a break in each curve between the regions where the quantization conditions 
for sufficiently small and sufficiently large values of m have been used. 





FIG. 4: Plots for the 2pa state of the ion of (a) C versus ri2 and (b) C versus ri2, when C and C are determined as 
functions of ri2 from the requirement that the first-order phase-integral results coincide with the numerically exact results 
obtained by Murai and Takatsu [13,14]. There is a break in each curve between the regions where the quantization conditions 
for sufficiently small and sufficiently large values of ri2 have been used. 



